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SUMMARY 

An algorithm has been developed for time-dependent forced convective diffusion-reaction 
having convection by a recirculating flow field within the drop that is hydrodynamically coupled at 
the interface with a convective external flow field that at infinity becomes a uniform free-streaming 
flow. The concentration field inside the droplet is likewise coupled with that outside by boundary 
conditions at the interface. A chemical reaction can take place either inside or outside the droplet, 
or reactions can take place in both phases. 

The algorithm has been implemented, and for comparison results are shown here for the 
case of no reaction in either phase and for the case of an external first order reaction, both for 
unsteady behaviour. For pure interphase mass transfer, concentration isocontours, local and average 
Sherwood numbers, and average droplet concentrations have been obtained as a function of the 
physical properties and external flow field. For mass transfer enhanced by an external reaction, in 
addition to the above forms of results, we present the enhancement factor, with the results now 
also depending upon the (dimensionless) rate of reaction. 

INTRODUCTION 

There are many industrial and environmental processes in which two-phase fluid-liquid 
systems are in use. Gases may be dispersed as bubbles in liquid phases, such as occurs in bubble 
columns and sparged vessels. Liquids may be dispersed in gases, such as occurs in scrubbers and 
in the atmosphere. And a liquid that is immiscible or partially miscible in another liquid may 
be dispersed in a liquid-liquid spray column extractor or reactor. The design of such systems 
may involve heat transfer, either intentionally or incidentally, but the widest range of applications 
involves mass transfer. 

Interphase mass transfer may proceed into or out of the dispersed phase. One (or more) 
chemical reaction(s) may take place in either the dispersed or the continuous phase in order to en- 
hance the rate of mass transfer. In two-phase reactions, certain of the reactants may be transferred 
from one phase into the other, where the reaction takes place, and the reaction products may then 
be transferred back into the first phase. Reactions may also occur in both phases. 

Because of the finite, generally small volume of each drop or bubble, interphase mass transfer 
unaccompanied by chemical reaction is inherently unsteady, regardless of the direction of mass 
transfer. Even if there is a reaction that admits of a steady state in the drop or bubble, unsteady 
behavior may nevertheless be of practical even primary importance. 

The continuous phase is inevitably in motion relative to the dispersed phase, and for clean 
systems (containing for instance no surface active agents) the motion in the two phases will be 
hydrodynamically coupled. 


291 



The general formulation could take into account either linear or nonlinear chemical reactions 
in either or both phases, and it could incorporate any velocity field that can be expressed as a 
function of r and i?. Neither are variable properties excluded. 

In order to demonstrate the utility of the algorithm, for concreteness we undertake the 
mathematical description of a liquid-liquid system in which a first order reaction takes place in 
the continuous phase. The dispersed phase is sufficiently dilute that the droplets which sediment 
(either falling under their weight or rising because of buoyancy) may be assumed isolated in an 
infinite medium, both with regard to fluid mechanics and to diffusion and reaction. The droplets 
are taken small enough that interfacial tension dominates shape effects and they are spherical. Al- 
though the approach we take and the methods we use do not require that that viscosity dominates 
flow effects and that the velocity fields have low Reynolds numbers, we nontheless consider the 
hydrodynamically coupled Hadamard - Rybczinsky profile for circulation within the droplet driven 
by an external velocity field that becomes a uniform streaming flow far from the droplet. Physi- 
cal and chemical properties are assumed constant, which would be the case for dilute isothermal 
systems, and we thus analyze interphase mass transfer for the forced convective diffusion-reaction 
single-drop system. We investigate specifically the role of the reaction rate, as measured by an 
appropriate Damkohler number, the solubility of the solute in the phases, as expressed by the lin- 
ear distribution coefficient (Henry’s law), the ratio of convection to diffusion, as measured by the 
Peclet number, and the ratio of the viscosities and that of molecular diffusivities of the two phases. 


GOVERNING EQUATIONS 


The dimensionless forced convective diffusion-reaction equations governing the solute con- 
centrations in the drop (0 < r < 1) and the continuous (1 < r < oo) phases, i = 1,2, respectively, 
can be represented in the form 


0cM 

dt 


+ K& vW.VeW = - jfWcW 


( 1 ) 


where i = 1 corresponds to the internal domain 0 < r < 1, and z = 2 to the external one 1 < r < oo. 

The dimensional parabolic partial differential equations have been rendered dimensionless 
using the droplet radius R as the characteristic length scale. The concentrations are measured in 
units of the initial driving force, 




(») gft) - Hc a 
Co R c Q 0 


; i = 1 , 2 , 


( 2 ) 


in which 




( 3 ) 


1, z = 1 
H, i = 2 , 

with H the Henry’s ”law” distribution coefficient, and 

t = 0 : c (1) = Co, c (2 > = Coo. (4) 

The characteristic time scale can be selected, for example, on the basis of the fastest physical 
or chemical process, occurring in the system, viz., 

r. = min(r£) nv , , r r ^ , i = 1,2), (5) 
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in which 


r« 

/ conv 


* t !& = — t« 

/WOOff-’ djff 2?W’ rxn 



( 6 ) 


The diffusivities and rate constants for the first order chemical reactions are denoted by !?(*) and 
JfeW, respectively, and the iT’s represent different combinations of standard dimensionless parameters 
for different choices of r,, as indicated in Table 1. 

Although our numerical implementation of the algorithm requires only that the velocity 
fields in the two phases be separable, we have selected the Hadamard • Rybczinsky solution for 
the convecting velocities in the dispersed and continuous phases to establish connections with 
earlier research [7, 8, 9]. In this instance, the characteristic velocity in each phase, with Uoo the 
freestreaming uniform flow at infinity, is taken as 


U? ) = fM(n)U 00) i = 1,2, (7) 

in which 

'"M ■ 20T Tr / <!, M = 1 • (•) 

with the viscosity ratio 


M = M ( V 2) - (9) 

The equations (10) are the ones used in the sequel, reflecting the selection of as the 
unit of time: 



. = 1 , 2 . 


with X = cost?, subject to the boundary conditions at the droplet interface, 


r = 1 : 


and at the limits of the overall domain, 



r = 0 : < oo 


( 11 ) 


( 12 ) 


r — * oo : c< 2 > —* 0 

Periodic boundary conditions in the angle variable, 


(13) 
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3c(‘> 


= 0, i = 1,2, 


(14) 


d=0,ir 

are satisfied automatically after introduction of the new independent variable A. 

The concentrations are subject to the initial conditions: 

t = 0: cW = 1, cW = 0 (15) 

The actual direction of mass transfer may be out of or into the drop, depending upon the 
driving force (co - i?Coo), even though the formulation of the problem suggests transfer from the 
droplet. 

The opposite direction of mass transfer in the actual problem would lead to the appearance 
of the inhomogeneous part in the reaction terms in (10) (but only when the corresponding ^ 0). 

THE ALGORITHM 

The problem is linear, and we use the Galerkin spectral method for the spatial discretization. 
The advantages of this method are well known [1, 2]. 

Boundary conditions at the origin of the droplet (12) and at infinity (13) and the symmetry 
boundary conditions (15) are implemented by the Lanczos tau-method [1, 2, 14]. 

We express the unknown functions c^(t, A,r) in a customary manner, 

c (i >(r,A,r)= £ cW(T,r) J>„(A), i = 1,2, (16) 

m=0 

in which the P m (A) are the Legendre polynomials of order m and the unknown coefficient functions 
Cm (t, r) are termed ’’radial functions” for brevity in the sequel. 

The discretization in the radial direction is performed in somewhat different ways for the 
internal and external domains. 

Using equation (10) for mass transfer inside the droplet ( i = 1), it is a simple matter to 
show that functions Cm\r, r) obey the following restrictions : 


c f (1) ( r » r=: °) = 0. 1 # 

(17) 

dc (1) 

» =°. 
r=0 

(18) 

c 2 ^(t, r) - even function of r 

k = 0,1... . 

2 ^ +1 (t, r) - odd function of r 

(19) 


On the basis of these restrictions, the radial functions inside the droplet were approximated 
by a series in even Chebyshev polynomials: 

c m ) ( r > r ) = 6 m , 0‘ q o(t) + rKm '^2 ^m,n( T )^2n-2(r), m = 0,l,...,M, (20) 

n= 1 

in which the T p (r) are Chebyshev polynomials of the first kind of order p, and 
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«2j = 2, = 


( 21 ) 


K\ — 1, fyj+i — 3, j — 1,2,... (22) 

Using (20) we automatically satisfy boundary condition (12) and avoid the singularity at 
the origin of the drop; the function <*o(t) represents the value of the concentration at the origin. 

Such an expansion on the interval 0 < r < 1 is valid as the even Chebyshev polynomials 
form a complete set for the type of functions considered [10]. 

The use of half the commonly used interval [-1,1] permits us to double the highest order 
of the polynomials used, leaving the number of terms in the series unaltered. 

The nonuniformity of the distribution of nodes in the spectral method (their number in 
close proximity to the surface is higher than near the origin) matches the physics of the problem 
as the concentration gradient near the interface is much larger. 

For the semi-infinite external domain we implement the widely used procedure of truncat- 
ing it at an appropriately large radius Too, far enough from the interface to make negligible the 
disturbance introduced by truncation. The boundary condition at infinity (13) is now imposed on 
this artificial boundary. It could be imposed as ’’hard”, ’’soft” [12] or ’’behavioral” [1, 13]. We use 
the ’’hard” one, 


r = Too ■ c< 2 ) = 0 (23) 

because it immediately results in the original boundary condition (13) if — ► oo. 

It is necessary to realize that by doing this we are changing the physical sense of the problem. 
The decrease to zero of the concentration infinitely far from its source is caused physically by the 
spreading of a finite amount of the species over an infinite spatial volume. After introduction of 
the boundary sphere at r = r^, we model this decrease by imposing what amounts to an infinitely 
fast heterogeneous reaction on the artificial boundary r^. The only justification for this is an a 
posteriori one, viz., by checking that the increase of r^ does not alter the solution in the vicinity 
of the droplet and in particular the interphase mass transfer. 

Our computations have confirmed this and show that when is chosen sufficiently large 
the choice of the particular type of boundary conditions mentioned above does not influence the 
resultant concentration distribution in regions where its value differs significantly from zero. 

The domain 1 < r < r^ is mapped onto the interval — 1 < z < 1 in such a way that the 
point z = 1 matches r = 1 and the point z — - 1 matches r = r „>• Among the wide variety of 
possible mappings two are used more often than others, the exponential and rational ones [1, 2]. 
A comparison by Grosch and Orszag [11] has shown that the latter mapping has some advantages 
over the former. 

Specifically, we use 


z = 


r — (1 + 6) 


I 1 -') ' 


(24) 


where S is the parameter representing the distance between the droplet surface and point mapped 
into z — 0. It is worth mentioning that we have also implemented the exponential mapping and 
could find no advantages for it over the rational mapping. 

The radial functions in the external domain are expanded as 
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N{i) 

C^ ) ( T » z ) = Yj ^mU T ) Zn(z), m = 0, 1, . . . , AT, (25) 

n= 1 

where the Z n (z), n = 1,2, are linear combinations of Chebyshev polynomials, each satis- 
fying the boundary condition following from (13): 

Z n (z = - 1) = 0, n=l,2,...,J\r( 2 >. (26) 

We take 


Zvkiz) - T 2 *(z) - 1 

(27) 

Zik-\{ z ) — T2k-i( z ) + 1 

Thus, we reduce the system of partial differential equations for two initially unknown func- 
tions c^\r, A, r) and cW(r, A, r) to a larger system of ordinary differential equations in r, for 

a o( r )> ^? ni . m ~ 0,1,..., M, (28) 

= n 2 = 1, 2, . . . , JV< 2 >. 

The total number of these unknown functions is 1 + ( M + 1) (N^ + N^). 

In order to obtain equations for these functions we use the conventional Petrov - Galerkin 
method, i.e., the basis functions are taken as the test functions [2]. We define two inner products: 



(29) 


(30) 


Forming by (29) the inner product of (10) for i = 1 with the test functions 


Po(\)T 0 (r) } P m (A)r' tm T 2ni _ 2 ( r), ms 0,1, 


,M, n 1 = l,2,...,W( 1 )-l, 


(31) 


and by (30) the inner product of (10) for t = 2 with the test functions 

PmWZ n3 {z), to = 0,1,..., AT, n 2 = 1,2,. .. , — 1, (32) 

we obtain two vector equations 



= (-At> + Kf - K& B<‘- r >) • b(‘> , i = 1,2. (33) 


Here A«, B^ c ), B^), B<‘- r ) are {1 + (M + 1)(JV<*) - 1), 1 + (AT + l)tf «} matrices, 
- {1 + ( M + 1)AT(*)} are the vectors of inhomogeneous terms, and <f>(*)(r) - {1 + ( M + 1 )JVW} 
are the unknown vectors, 


<£ (1) ( r ) = ( «o» 4> ( o,l 


■<t> 


,(i) 

o,Ar(o> 






(i) 

Jlf.Ah 1 ) 


(34) 
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(35) 


*<’>(■ 




■ ■ > ^M,NW ) T - 


The remaining 2 (M + 1) equations are derived from the boundary conditions (11). 

Upon substituting (20) and (25) into (11), multiplying by P m (A), m = 0,1 ,..., M and 
integrating A from -1 to 1, we obtain two sets of M+l linear algebraic equations: 




(36) 


H-D-SW-t^ = sW-<t>( 2) , (37) 

where QM, S<‘) are {(M + 1), (Af + 1)(1 + N (*))} matrices, i = 1,2. 

By expressing and <f>^ N w , to = 0,1,..., AT, using the system (36)-(37) and 

substituting in the system (33), we arrive finally at the system of 1 + ( M + l)(JV( l ) + N& - 2) 
linear ODEs: 


A^- = (B (c) + B^ + B( r ))-d> + b. (38) 

dr 

The constant matrices B^, B ^ and B^ r ) correspond respectively to the convective, 

diffusive, and reactive terms in the original equation (10), b is an {1 + (M + 1 )(N^ + - 2)} 

constant vector, and 4> (r) is the vector of unknown functions 


<j) = (ao, ^>0,11 • • • ^o,N( l )-i > ^o,i > • • • ^o,Af( J )-i ’ • • ' » 


1» 




kW 


and not simply a concatenation of vectors <f> ^ and <j> 

The matrices A, B W and B( r ) are block-diagonal. They all have M+l nonzero square 
{jy(i)-|-jV'( 2 )-2 , N^+N^-2)} matrices on their main diagonals and their first l+(N^+N^-2) 
elements in the first row and the first column are nonzero. 

The matrices B( c ) that result from transforming the convective terms also have block struc- 
ture with the same block sizes. However, they are no longer block-diagonal and the amount of 
nonzero block-diagonals depends on the velocity fields v(‘), i = 1,2. The higher the degree of A 
that is involved in the velocity field expressions, the greater will be the coupling between the radial 
functions of different orders. And the increase of the order of this coupling leads to a corresponding 
increase in the number of nonzero block diagonals in B( c ). 

For the Hadamard - Rybczinsky field, for example, these matrices are block-tridiagonal, but 
for the velocity field in [15], valid for higher Reynolds numbers, it would be block-pentadiagonal. 

The discontinuous initial conditions (15) are not appropriate for computations. Instead, 
we used the analytical solution for the pure diffusion case (no convection, no chemical reaction) 
derived in [16]. The concentration distributions for very small time values were expanded over our 
basis functions 72 „- 2 ( r ) and Z n , n = 1,2,... to initialize the computations, and the coefficients 
obtained were used as initial conditions for Qo( r ), $m?n, , and ^>m,n 2 , to = 0, 1,.. ., M, ni = 
1,2, . . ., jyf 1 ) — 1, n a = l,2,...,J\r< a >-l. 

For time discretization of the system (38) we used the first-order backward Euler method. 
Defining 4> n as vector <t> at the n-th time step of magnitude At and 


B = B (c) + B (d) + B (r) , 

system (38) can be rewritten as 


(40) 


297 



where 


(A - ArB)-A«£ n+1 = ArB-<f> n + Arb , 


(41) 


A* n+1 = </> n+1 - <f> n . (42) 

Every time step system of linear equations (41) was solved by regular Gauss elimination 
(preceded by LU decomposition) with the following iterative refinement [3]. The matrix on the 
1 left side of (41) has the same structure as the matrix B; as mentioned, it is block-tridiagonal 
for the Hadamard-Rybczinsky velocity field. Our attempts to apply block-elimination methods (in 
particular, the block Thomas algorithm [4]) failed, presumably because block LU factorization does 
not involve pivoting, which is essential when diagonal dominance does not occur (which is the case 
for high Peclet numbers). 

We considered the matrix on the left side of (41) as a banded one with bandwidth 1 + 
3 (WW + tf( 2 ) - 2). 

As long as this matrix depends on the time step and its factorization is a time-consuming 
process, only two values of the time step were used for each run. A smaller one was used for an 
initial time period and an another one for the subsequent time range. 

The numbers of terms in series (16), (20), and (25) depend on the steepness of the con- 
centration gradients and were different for different values of Peclet and Damkohler numbers. The 
maximum numbers used were M = 87, N = 25, = 97. 

As is well known [1, 2], an increase in the number of terms in a spectral series (especially 
in the series in Chebyshev polynomials) leads to very high condition numbers for the resulting 
system of linear equations. This was alleviated by using double precision in all computations and, 
as mentioned above by application of the iterative refinement to the solution obtained with the 
Gauss elimination procedure. 


QUANTITIES OF INTEREST 

The most practically interesting quantity in extraction problems is the amount of material ex- 
tracted by a particular instant in time. For the problem under consideration (i.e., when species are 
extracted from the droplet) this can be conveniently characterized by the time-dependent average 
dimensionless concentration of species remaining in the drop: 

= J/cM(r,X,r) d Xdr ( 43 ) 

This quantity changes in time as a result of mass transfer out of the droplet. The local and 
surface average rates of this transfer are characterized by corresponding mass transfer coefficients, 
the quantities which when multiplied by the driving force give the respective mass transfer rate. 
The nondimensional mass transfer coefficient is usually referred to as the Sherwood number Sh, 
which is analogous to the Nusselt number in heat transfer problems. 

Different kinds of Sherwood number can be introduced, depending on the driving force on 
which they are based and the domain to which they are related. 

For the problem of single-drop extraction, the instantaneous driving force F'( dr ) for mass 
transfer is the difference between the concentration of the transferring species in the droplet and 
that far away from it, taking into account the step change of the concentration at the interface due 
to solubility, 
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(44) 


F (dr) = |(1> _ H - Coo t 

where is the dimensional average concentration of species in the droplet. 

Often the Sherwood number is based on the maximum possible (in our case, initial) driving 

force: 

= co - Hc^ . (45) 

Here we consider only the external Sherwood number, i.e., the nondimensional rate of 
transfer of species from the external side of droplet surface into the external flow. 

The local and average external Sherwood numbers defined on the basis of the maximum 
driving force are respectively: 


Sh\ OCt o = —2 H'D 


dcM 

dr 


r = 1 


and 


Sh 0 


- b - b L 


1 dcM 


dr 


d\. 


r=l 


Corresponding values based on the instantaneous driving force are: 

Sh ioc,o 


Sh\ oc — 


Sh = 


d 1 ) 

Shp 
7*0 • 


(46) 


(47) 


(48) 


(49) 


The che mi cal reaction in the external region increases the rate of the extraction, and this 
increase is characterized by the enhancement factor, which is the ratio of the corresponding mass 
transfer rates [5]: 


E ShjD^ £ 0 ) 

Sh(Da { fl = 0) 


(50) 


COMPUTATIONAL RESULTS AND DISCUSSION 

The results of the computations to be presented cover the following ranges of parameters: 

0.25 < D < 4 , 

0 < Pe (2) < 500 , 

0 < Dafi < 1000 , 

H = fi = 1 . 

The characteristic time scale was chosen as 

r -r (2) 

T * - r diff > 


(51) 
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which is just the Fourier number based on the diffusion coefficient of the external fluid. The times 
appearing on the plots are expressed in these units. The values of Peclet number Pe and Damkohler 
number Da presented on the plots correspond to Pe and Da$ , respectively. 

To illustrate qualitively the process of pure mass transfer (no reaction) from the droplet, 
we present in Figures 1-3 the isocontours of constant species concentration at different times for 
various levels of external convection (Pe^ = 10, 200, 500, respectively). The well known and 
intuitively expected increase of mass transfer with increasing convection is apparent. 

The influence of internal circulation on the development of the mass transfer process is 
illustrated in Figures 4-5 where we present the isoconcentration contours for the same external 
Peclet number (Pe< 2 ) = 500) and different ratios of internal and external diffusivities ( D = 0.25 
and D = 4.0). 

For D = 0.25 the internal convection is much stronger, in the sense that the value of Pe^ 1 ) 
is larger. As a consequence the concentration isocontours inside the droplet lie close to the internal 
streamlines, a result already obtained numerically by Johns and Beckmann [7] for the special case 
of mass transfer resistance solely inside the droplet. The coincidence of internal isocontours with 
internal streamlines also constituted the basic assumption of Kronig and Brink’s model of mass 
transfer in a circulating drop [6]. From a simple comparison of the isocontour levels in Figures 4 
and 5 alone one infers that the mass transfer from a droplet for D=4 is much more intense than 
for D=0.25. The reason that the internal Peclet number Pe is greater for D = 0.25 is not that 
the internal circulation is greater, for it is not (fi = 1), but that the internal diffusivity is smaller. 
Nonetheless, it is customary for brevity to describe an increase in Peclet number as an increase 
in convection, rather than the more lengthy but more accurate increase of the ratio of convection 
to diffusion. In this usage, one may phrase the conclusion drawn from Figures 1-5 as follows: 
increased convection outside the droplet increases the rate of extraction, but increased convection 
inside suppresses mass transfer. 

The influence of the external reaction rate on the concentration distribution is shown in 
Figures 2, 6 and 7. As could have easily been anticipated, an increase in Da$ results in faster 
extraction and an almost immediate disappearence of extracted species outside the droplet (almost 
no species present for Da$ = 100 in Figure 7). 

Figures 8 and 9 show the effect of reaction rate on the local Sherwood number. The values of 
Sh\ oc o go to zero with time for all values of the angle variable i?, although the surface distribution of 
Sherwood number based on the instantaneous driving force approaches a nonvanishing asymptote. 
An increase in the reaction rate thus results in a general increase of mass transfer and of values of 
the Sherwood numbers, but the temporal variation of values of local Sherwood numbers at different 
locations is less transparent, warranting further investigation. 

Figure 10 reflects the behavior of average Sherwood number Sh in time for different values 
of the external Peclet number Pe ( 2 ) for the no- reaction case. The oscillations of Sh were obtained 
computationally by different investigators including Johns and Beckmann [7] and Oliver and Chung 
[9], who solved the conjugate unsteady heat transfer problem, which is mathematically equivalent to 
the mass transfer problem under consideration when there is no chemical reaction involved. These 
oscillations are caused by the internal circulation, with the most detailed physical explanation being 
given by Brignell in [8]. Consequently, the period of these oscillations is smaller and the amplitute 
greater the higher the Peclet number is. The stronger convection also leads to a higher mass transfer 
rate as it creates the thinner diffusion boundary layers on the both sides of the droplet surface. 

Figure 11 illustrates the influence of the rate of external chemical reaction on the average 
Sherwood number. The plots confirm the conclusions made above on the increase of the rate of 
extraction with the increasing external convection and rate of external chemical reaction. 
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In a more apparent way, this is reflected in Figure 12, where the decrease of the average 
droplet concentration with time is shown. From this picture we can also deduce the very important 
conclusion that an increase in the reaction rate beyond some specific value will not benefit the 
extraction results (the diflfences between the average droplet concentration for Da^jj values of 300 
and 1000 are quite small). 

Figure 13 shows the effect of reaction rates on the values and temporal evolution of the 
enhancement factor E. The oscillations here are the consequences of internal circulation, the same 
as for the corresponding average Sherwood number on Figure 10. The values of E corresponding 
to the same reaction rate are higher for smaller Pe (lower convection). A possible explanation 
could be that the corresponding values of Daff\ which are just the ratios of Da$ and Pe^ 2 \ are 
smaller for higher Pe^ 2 \ 

In closing, we want to emphasize that the purpose of this article has been to present the 
numerical algorithm we have developed and to illustrate some of the kinds of results that can be 
obtained for this concrete situation. Our subsequent articles will include additional results for other 
classes of forced convective diffusion-reaction problems for single drops, as well as more detailed 
analyses of these and other results. 
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NOMENCLATURE 




Co 

Coo 

c(‘) 
c<0 
I>( 0 
D 

Pag 

E 

/(0 

A« 

M 

Af(‘) 

Pc(‘) 

r 

R 

£f<‘> 

t/oo 

v(*') 


1? 

A 

M 1 
r 

Teonv 


(«) 


T (<) 

T diff 

Trxn 


•dimensional value of the uniform concentration in the origin of the droplet at t = 0 

-dimensional value of the concentration far from the drop 

-dimensional concentration in the t-th domain, t = 1,2 

-dimensionless concentration in the t-th domain, t = 1,2 

-molecular diffusivity of the solute in the fluid in the t-th domain, t = 1,2 

-molecular diffusivity ratio, D^/D^ 


-first Damkohler number in the t-th domain , 


OR 

ls!±Bi i 


•second Damkohler number in the t-th domain, K A 

D' 1 


t = 1,2 
= 1,2 


-enhancement factor, eq.(50) 

-factor showing the viscosity ratio dependence of the 
velocity scale in the t-th domain, t = 1,2 
-distribution coefficient (Henry ’’law” constant) 

•chemical reaction rate constant in the t-th domain, t = 1,2 
-highest order of the Legendre polynomials used in the expansion 
in the angular direction 

-number of terms in the expansion of radial functions in the t-th domain , t = 1 


,2 


-Peclet number in the t-th domain, 


2UoofM(n)R 

i)(‘) 


t = 1,2 


-dimensionless radial coordinate 
•droplet radius 
•dimensional time 


-characteristic velocity scale in the t-th domain, t = 1, 2 
-velocity of the flow at the infinity 
•velocity field in the t-th domain nondimensionalized 
by the corresponding velocity scale U»\ t = 1,2 
-polar angle in spherical coordinate system 

= COST? 

-molecular viscosities ratio, /t^ 1 )/^ 2 ) 

•molecular viscosity of the fluid in the t-th domain, t = 1,2 
-dimensionless time 


-convection time scale in the t-th domain, 


o 


-, i=l,2 


-diffusion time scale in the t-th domain, t = 1,2 

-chemical reaction time scale in the t-th domain, ^y, t = 1,2 
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Fig. 9 Distribution of the values of local Sherwood number Sh at Pe=200 for different values of Da 










Fig. 1 1 Effect of reaction rate on average Sherwood number Sh for Pe=200 and Pe=500 
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Fig. 12 Time evolution of average droplet concentration 
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